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A  free  vibration  analysis  of  a  polymer  electrolyte  membrane  fuel  cell  (PEMFC)  is  performed  by  modelling 
the  PEMFC  as  a  20  cm  x  20  cm  composite  plate  structure.  The  membrane,  gas  diffusion  electrodes,  and 
bi-polar  plates  are  modelled  as  composite  material  plies.  Energy  equations  are  derived  based  on  Mindlin’s 
plate  theory,  and  natural  frequencies  and  mode  shapes  of  the  PEMFC  are  calculated  using  finite  element 
modelling.  A  parametric  study  is  conducted  to  investigate  how  the  natural  frequency  varies  as  a  function 
of  thickness,  Young’s  modulus,  and  density  for  each  component  layer.  It  is  observed  that  increasing  the 
thickness  of  the  bi-polar  plates  has  the  most  significant  effect  on  the  lowest  natural  frequency,  with  a 
25%  increase  in  thickness  resulting  in  a  17%  increase  in  the  natural  frequency.  The  mode  shapes  of  the 
PEMFC  provide  insight  into  the  maximum  displacement  exhibited  as  well  as  the  stresses  experienced 
by  the  single  cell  under  vibration  conditions  that  should  be  considered  for  transportation  and  stationary 
applications.  This  work  provides  insight  into  how  the  natural  frequencies  of  the  PEMFC  should  be  tuned 
to  avoid  high  amplitude  oscillations  by  modifying  the  material  and  geometric  properties  of  individual 
components. 
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1.  Introduction 

Polymer  electrolyte  membrane  fuel  cells  (PEMFCs)  may  be  sub¬ 
ject  to  vibrations  under  dynamic  situations  found  in  transportation 
applications  as  well  as  stationary  applications  near  heavy  traffic 
and  rail  transport.  Passenger  vehicles  generally  experience  vibra¬ 
tions  in  the  range  of  8-1 6  Hz  due  to  the  unevenness  of  the  road 
and  the  oscillation  of  the  axel  and  wheel  with  the  suspension  sys¬ 
tem  [1].  PEMFCs  may  be  employed  to  power  auxiliary  devices  in 
semi-trailers,  which  typically  experience  vibrations  in  the  range  of 
0.9-5.8Hz  during  highway  driving  conditions  [2].  Buildings  near 
busy  roads  are  also  subject  to  vibrations  due  to  nearby  traffic, 
which  cause  vibrations  between  5  and  25  Hz  [3].  When  placed  in 
applications  such  as  these,  the  PEMFC  may  vibrate  at  an  excita¬ 
tion  frequency  within  the  band  of  its  natural  frequencies  and  cause 
the  PEMFC  to  vibrate  in  resonance  with  high  amplitude.  Therefore, 
it  is  critical  to  identify  the  natural  frequencies  that  a  PEMFC  may 
experience  in  the  context  of  the  excitation  frequencies  expected 
in  these  applications.  Vibrating  at  resonance  frequency  can  lead 
to  the  initiation  and  acceleration  of  defect  formation,  which  may 
ultimately  result  in  operational  failure.  Vibrations  may  exacerbate 
defects  such  as  pinholes,  cracks,  and  delamination,  which  can  result 
in  fuel  crossover,  performance  degradation,  and  reduced  durabil¬ 
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ity  [4,5].  Vibration  characteristics  are  required  to  understand  the 
vibration  behaviour  of  PEMFC  components  such  as  the  membrane, 
catalyst  layer,  gas  diffusion  layers,  and  bi-polar  plates. 

Vibration  tests  on  PEMFC  stacks  and  mechanical  non-linear 
behaviour  modelling  with  neural  networks  and  finite  element  mod¬ 
elling  has  been  reported  by  several  authors  [6-13].  Rouss  et  al. 
[7]  test  a  fuel  cell  stack  under  vibrations  to  simulate  the  opera¬ 
tional  conditions  of  an  aircraft.  They  characterize  the  mechanical 
behaviour  of  the  stack  by  determining  the  natural  frequencies,  and 
they  find  that  the  PEMFC  behaves  as  a  multi-body  non-linear  sys¬ 
tem.  Each  component  of  the  fuel  cell  exhibits  distinct  frequency 
responses,  in  contrast  to  behaving  like  a  rigid  body.  In  a  sepa¬ 
rate  work,  the  authors  develop  a  three-dimensional,  multi-input 
and  multi-output  artificial  neural  network  (ANN)  model  based  on 
experimental  findings  from  PEMFC  vibration  tests  [8].  They  pro¬ 
pose  a  fuel  cell  neural  network  model,  which  can  be  used  for 
monitoring  purposes  to  detect  abnormalities  in  the  mechanical 
behaviour  of  a  similar  fuel  cell  stack  placed  under  vibrating  con¬ 
ditions. 

The  mechanical  non-linear  response  of  a  PEMFC  is  characterized 
using  acceleration  data  from  vibration  experiments  by  Rouss  et  al. 
[  1 0].  They  place  accelerometers  on  the  bi-polar  plate  with  and  with¬ 
out  structural  supports  and  find  similarities  between  the  results 
from  a  simple  single  degree  of  freedom  (SDOF)  system  and  from  a 
complex  multi-physics  system  based  on  a  library  of  non-linear  sig¬ 
natures.  From  temporal  signals  measured  with  one-dimensional 
and  three-dimensional  accelerometers,  they  find  that  the  PEMFC 
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Nomenclature 

A  cross-section  area 

cii  coefficients  of  mode  shape 

Db  the  material  property  matrix  for  bending 

Ds  the  material  property  matrix  for  shear 

Dy  the  coefficients  of  the  bending  and  shear  stiffness 
matrices 

Ek  modulus  of  elasticity  of  the  respective  layer 

h  thickness  of  the  PEMFC 

L  length  of  the  PEMFC 

[I<]  assembled  stiffness  matrix 

I<e  element  stiffness  matrix 

[M]  assembled  mass  matrix 

Me  element  mass  matrix 

n  the  number  of  nodes  per  element 

nc  the  number  of  component  layers 

Nj  the  shape  functions 

t  time  variable 

T  the  kinetic  energy 

tk  thickness  of  the  respective  layer 

U  the  strain  energy 

V  the  volume  of  the  PEMFC 

w  displacement  in  the  thickness  direction 

{w}  displacement  vectors 

{w}  acceleration  vectors 

x  longitudinal-direction  of  the  PEMFC 

y  width-direction  of  the  PEMFC 

z  thickness-direction  of  the  PEMFC 

crb  the  bending  stress  tensor 

cr $  the  shear  stress  tensor 

£b  the  bending  strain  tensor 

ss  the  shear  strain  tensor 

0X  rotation  about  x-axis 

0y  rotation  about  y-axis 

pk  density  of  the  respective  layer 

a  shear  correction  factor 

X  Eigen  value 

c on  natural  frequency 

vk  Poisson’s  ratio  of  the  respective  layer 

§, rj  local  coordinates 


exhibits  dry  friction  behaviour,  which  they  attribute  to  the  con¬ 
tact  between  the  bi-polar  plates  and  tie  rods.  To  investigate  the 
dynamic  behaviour  of  a  non-linear  mechanical  system,  Rouss  and 
Charon  [11]  measure  the  response  of  a  non-linear  mechanical  sys¬ 
tem  subjected  to  random,  shock,  and  swept  sine  excitations  on  a 
vibrating  platform.  They  employ  these  results  to  validate  and  train 
their  numerical  model. 

Rajalakshmi  et  al.  [6]  subject  a  PEMFC  stack  to  vibrations,  which 
include  random  and  swept-sine  excitations  on  a  vibrating  platform 
in  three  axes.  Although  changes  to  the  mechanical  integrity  of  the 
stack  are  not  detected,  they  find  a  compression  force  release  at  the 
bolts.  Ozgur  [14]  performs  modal  analysis  with  experimental  vibra¬ 
tion  investigations  on  an  automotive  fuel  cell  stack  module  in  order 
to  design  a  mounting  bracket.  They  employ  an  electro-dynamic 
shaker  test  apparatus  to  identify  the  appropriate  mounting  system 
under  resonance  frequency  considerations. 

Betournay  et  al.  [12]  experimentally  investigate  the  effects  of 
mining  conditions  on  the  performance  of  a  PEMFC.  They  determine 
the  effects  of  shocks  and  vibrations  on  the  performance  of  PEMFC 
and  the  effects  of  mineral  and  diesel  particulate  matter  on  the  phys¬ 
ical  reliability/integrity  of  the  fuel  cell  stack.  They  also  evaluate  the 
physical  reliability/integrity  of  the  fuel  cell  stack  when  mounted 


over  the  rear  wheel  of  a  mine  loader  chassis,  which  applied  shocks 
and  vibrations  to  the  fuel  cell. 

Privette  et  al.  [15]  conduct  vibration  studies  to  qualify  a  PEMFC 
against  US  military  vibration  and  shock  standards.  PEMFC  stack 
testing  has  been  recently  studied  in  the  framework  of  a  European 
Union  harmonized  fuel  cell  testing  protocol  [16],  which  is  com¬ 
prised  of  55  European  partners.  Fuel  cell  testing  protocols  include 
the  application  of  vibrations  and  shocks  with  6  degrees  of  freedom 
at  a  frequency  of  up  to  250  Hz. 

To  the  authors’  best  knowledge,  free  vibration  analysis  using 
the  finite  element  method  has  yet  to  be  performed.  In  this  work, 
the  finite  element  model  for  simulating  the  PEMFC  as  a  composite 
material  is  presented.  Free  vibrations  are  simulated  for  the  compos¬ 
ite,  which  consists  of  a  membrane,  gas  diffusion  electrodes  (GDEs), 
and  bi-polar  plates.  A  parametric  study  is  conducted  to  investigate 
how  the  natural  frequency  varies  as  a  function  of  thickness,  Young’s 
modulus,  and  density  for  each  component  layer. 

2.  Finite  element  formulations 

In  this  work,  the  finite  element  method  is  employed  for  the 
free  vibration  analysis  of  a  PEMFC,  with  degrees  of  freedom  for 
the  transverse  displacement  and  rotation  with  respect  to  the  x- 
and  y-directions.  The  PEMFC  is  modelled  as  a  symmetrically  lam¬ 
inated  composite  plate,  composed  of  a  membrane,  gas  diffusion 
electrodes  (GDEs)  and  bi-polar  plates,  which  are  considered  plies 
of  the  composite  plate.  The  GDE  is  the  gas  diffusion  layer  combined 
with  the  catalyst  layer.  The  bi-polar  plate  is  modelled  as  a  solid 
graphite  layer,  without  flow  fields.  Since  the  PEMFC  is  modelled 
as  a  composite,  interfaces  between  the  layers  are  not  considered 
here.  For  the  free  vibration  analysis  of  a  PEMFC,  clamped  boundary 
conditions  are  applied  to  all  edges. 

Isoparametric  shape  functions  are  employed  for  the  plate  ele¬ 
ment  formulations.  The  generalized  transverse  displacement  and 
rotations  are  interpolated  as  follows  [17]: 


n 


w  =  r,)W, 

1=1 

(la) 

n 

ex  =  J2Ni&  iV*, 

i=l 

(lb) 

n 

9y  =  ri0y> 

J=1 

(lc) 

where  n  is  the  number  of  nodes  per  element.  represents 

the  shape  functions  of  a  bilinear  four-node  quadrilateral  element 
in  natural  coordinates,  based  on  a  local  coordinate  system  with 
independent  variables^  and  rj.  w,  0X  and  0y  are  the  transverse  dis¬ 
placements  in  the  z-direction,  and  rotations  of  the  normal  to  the 
mid-plane  with  respect  to  the  x  and  y-direction,  respectively.  wit 
0Xi ,  and  0yi  are  the  nodal  transverse  displacements  in  the  z-direction 
and  the  nodal  rotation  with  respect  to  x-direction  and  y-direction, 
respectively.  The  shape  functions  of  a  4-node  bilinear  quadrilateral 
element,  which  are  used  in  the  strain-displacement  relation  of  the 
energy  formulation  section,  are  as  follows: 

(2a) 

N2&»7)=  J(i  +£Xi-»?) 

(2b) 

n3(i,  to  =  J(i  +m+v) 

(2c) 

n4(I,»?)=  J(WXi  +  »?) 

(2d) 
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Mindlin’s  plate  theory  [  1 8  ]  is  utilized  for  the  energy  formulation, 
which  incorporates  the  transverse  shear  deformation.  The  funda¬ 
mental  assumption  is  that  the  plate  cross-section  remains  plane 
and  normal  to  the  deformed  longitudinal  axis  before  bending,  but 
does  not  necessarily  remain  normal  to  the  deformed  middle  plane 
of  the  plate  after  bending.  Taking  into  account  the  basic  assump¬ 
tions  of  Mindlin’s  plate  theory  for  a  PEMFC,  the  strain  energy  is 
written  as: 

u=l  fcrbsbdv+^  [asSsdV  (3) 

Jv  J  V 

where  crb  denotes  the  bending  stress  field  tensor,  sb  denotes  the 
bending  strain  field  tensor,  crs  denotes  the  transverse  shear  stress 
field  tensor,  and  es  denotes  the  transverse  shear  strain  field  tensor. 
V  is  the  volume  of  the  plate,  and  a  is  the  shear  correction  factor. 
Subscripts  b  and  s  denote  bending  and  shear,  respectively. 

The  linear  elastic  stress-strain  relations  in  bending  are  defined 
for  an  isotropic  material  as  the  following: 


(7b  —  Db£b 


(4) 


The  total  functional  energy  in  the  absence  of  work  from  external 
forces  is  expressed  as: 


n  =  U-T  (12) 

By  imposing  the  stationary  condition  ( 3/7/ 9w  =  0)  to  Eq.  (12), 
the  element  stiffness  matrix,  Ke,  and  mass  matrix,  Me,  are  calculated 
as  follows: 

I<e=  I  BTbDbkBbdA  +  a  I  BTsDskBsdA  (13) 

JA  Ja 


r  f™ 

tk 

0 

0  ■ 

Me=  /  w  y^Pk 

0 

tfc3/12 

0 

V  k=l 

0 

0 

t3k/u  \ 

NdA 


(14) 


By  substituting  the  local  coordinates  in  the  element  area, 
the  element  stiffness  and  mass  matrices  are  computed  with  the 
Gauss-Legendre  quadrature  numerical  integration  as  follows  [17]: 


I<e  = 


BThDbB„U  ]d^drj+a 


BTsDsBs\J]d^dr] 


(15) 


The  material  property  matrix  for  bending  of  a  symmetric  lami¬ 
nated  composite  plate  is  given  by  [19]: 


Db 


Du  £>12  £>16 

Di2  D2  2  D2  6 
Die  D26  Dee 


where 


D22  =  Du 


(5) 


(6a) 

(6b) 

(6c) 


Die  =  D26  =  0 


D66  =Yj 


nc  VkEk 
k=i  2(1  -V2k) 


(6d) 

(6e) 


Ek  denotes  Young’s  modulus,  vk  denotes  Poisson’s  ratio,  tk  denotes 
the  thickness  of  the  component  layers,  and  the  subscript  k  denotes 
the  respective  component  layer. 

The  linear  elastic  stress-strain  relationship  for  shear  is  defined 
for  an  isotropic  material  as: 


&s  —  Ds£s 


(7) 


The  material  property  matrix  for  the  shear  of  a  symmetric  lam¬ 
inated  composite  plate  is  given  by: 


Dc  = 


D77  0 

0  Dg  8 


where 

E/il 

T 

k=  1  1 


Ek 


~tk 


:=1  1  -  Vk 

The  maximum  strain  energy  of  the  plate  is  written  as: 


U=2 


l  f  [£b]T[Db][sb]dV+^  f  [es]r 
Jv  Jv 


[Ds][ss]dV 


The  kinetic  energy  for  a  plate  is  given  by: 


T  = 


if 


dA 


(8) 


(9) 


(10) 


(ID 


where  p  is  the  density,  h  is  the  thickness  of  the  plate,  and  A  the 
cross-sectional  area. 


N\j\d^drj  (16) 


where  |j|  is  the  determinant  of  the  Jacobian  matrix.  After  calculat- 
ing  element  stiffness  and  mass  matrices  using  a  numerical  program 
written  in  Matlab  [20],  the  matrices  are  assembled  together  to  form 
the  overall  stiffness  and  mass  matrices  for  a  complete  single  PEMFC. 
These  assembled  stiffness  and  mass  matrices  are  used  in  the  equa¬ 
tion  of  motion,  which  is  discussed  in  the  section  below. 


3.  Solution 


The  equation  of  motion  for  a  free  vibration  of  a  structure  is 
written  as: 

[M]{w}  +  [I<]{w}  =  0  (17) 

where  [M],  [JC],  {w}  are  the  mass  matrix,  stiffness  matrix,  and  dis¬ 
placement  matrices,  respectively.  To  find  the  natural  motion  of  a 
structure,  the  form  of  response  or  solution  is  assumed  as: 

{w(t)|  =  {Z}eia)t  (18) 

where  {Z}  is  the  mode  shape  (Eigen  vector),  and  co  is  the  natu¬ 
ral  frequency.  The  general  solution  is  a  linear  combination  of  each 
mode: 

{w(t)|  =  aifZile1'^  +  a2{Z2Wco  ^  +  a3{Z3}eico^  + . . .  (19) 

where  each  constant  (a,-)  can  be  evaluated  from  the  initial  condi¬ 
tions.  Substituting  Eq.  (18)  into  Eq.  (17)  yields: 

([K]-  <w2[M]){Z)ei®t  =  0  (20) 

The  above  equation  has  a  nontrivial  solution  if  |[K]  -  <w2[M]| 
becomes  singular.  In  other  words,  there  exists  n  number  of  solu¬ 
tions,  {c o\ ,  co%, . .  .co%),  which  satisfy  the  following  equation: 

\[K]-X[M]\[Z)  ={0}  (21) 

where  A  =  co2  is  the  Eigen  value  of  the  system. 

From  a  grid  sensitivity  study,  a  10  x  10  element  mesh  is  found 
to  provide  reasonable  convergence  and  is  utilized  for  all  simula¬ 
tions  presented  in  this  work.  The  lowest  three  natural  frequencies 
are  calculated  and  presented  in  Table  1.  Based  on  this  convergence 
study,  doubling  the  mesh  density  would  only  result  in  a  0.2%  dif¬ 
ference  in  the  calculated  first  natural  frequencies.  Therefore,  the 
mesh  density  employed  in  this  work  is  sufficient.  A  schematic  of 
the  discretized  PEMFC  is  shown  in  Fig.  1. 
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Table  1 

Comparison  of  natural  frequencies  (Hz)  of  a  PEMFC  for  different  element  meshes. 


Mode  number 

5  x  5-mesh 

6  x  6-mesh 

8  x  8-mesh 

10  x  10-mesh 

12  x  12-mesh 

15  x  15-mesh 

First  mode 

2328.80 

2318.10 

2310.70 

2308.54 

2307.73 

2307.36 

Second  mode 

4430.50 

4386.80 

4356.80 

4450.88 

4340.08 

4344.08 

Third  mode 

6064.00 

6009.20 

5966.70 

5965.83 

5956.67 

5954.95 

Fig.  1.  A  schematic  of  the  discretized  mesh  for  the  PEMFC  (a)  top  view  and  (b) 
three-dimensional  view. 


4.  Parametric  study 

A  parametric  study  is  performed  to  investigate  the  mechanical 
and  geometrical  property  effects  on  the  vibration  characteristics 
of  a  PEMFC  The  thickness,  Young’s  modulus,  and  density  for  each 
component  layer  are  increased  independently  in  increments  of  5% 
to  a  maximum  of  25%  of  the  base  case  value,  and  the  impact  on  the 
natural  frequency  of  the  first  mode  is  determined.  Mechanical  and 
geometrical  properties  for  the  base  case  simulation  are  obtained 
from  [21  ]  and  are  listed  in  Table  2. 

The  mode  shapes  for  this  base  case  are  shown  in  Fig.  2,  which 
provides  insight  into  the  maximum  displacement  exhibited  as  well 
as  the  stresses  experienced  by  the  material. 


4.1.  Influence  of  material  thickness 

Fig.  3  illustrates  the  impact  of  independently  increasing  the 
thickness  of  the  PEMFC  components  on  the  lowest  natural  fre¬ 
quency.  5%  increments  of  the  thickness  in  the  bi-polar  plate,  GDE, 
and  membrane  result  in  increases  of  thickness  on  the  order  of  1 0_1 , 
10-2  and  10-3,  respectively.  The  dominating  effect  of  the  bipolar 
plate  thickness  is  exacerbated  by  the  t3  term  in  the  material  prop¬ 
erty  matrix  (Eqs.  (6a-e),  (9)  and  (14)).  Therefore,  increasing  the 
thickness  of  the  GDEs  and  membrane  result  in  negligible  changes 
to  the  stiffness  and  mass  matrices  and  negligible  increases  to  the 
lowest  natural  frequency  (less  than  1%).  However,  increases  to  the 
thickness  of  the  bi-polar  plates  result  in  a  significant  linear  increase 
in  the  lowest  natural  frequency,  with  an  increase  of  approximately 
17%  when  the  thickness  of  the  plates  increases  by  25%. 


(C) 


0  0 


Fig.  2.  Normal  mode  shapes  for  a  PEMFC  considered  as  a  Mindlin’s  composite  plate 
(a)  first  mode,  (b)  second  mode,  and  (c)  third  mode. 


Fig.  3.  Effect  of  independently  varying  component  thickness  on  the  first  natural 
frequency. 
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Table  2 

Mechanical  and  physical  properties  of  the  PEMFC  for  the  base  case  simulation  [21  ]. 


Young’s  modulus  (GPa) 

Poisson’s  ratio 

Density  (kg  m-3) 

Thickness  (mm) 

Length  (mm) 

Width  (mm) 

Bi-polar  plate 

10 

0.25 

1800 

12.8 

200 

200 

GDE 

10 

0.25 

400 

0.28 

200 

200 

Membrane 

0.197 

0.25 

2000 

0.05 

200 

200 

Fig.  4.  Effect  of  independently  varying  Young’s  modulus  on  the  first  natural  fre¬ 
quency. 

4.2.  Influence  of  Young's  modulus 

The  effect  of  independently  increasing  the  Young’s  modulus  for 
the  bi-polar  plates,  GDEs  and  membrane  on  the  lowest  natural  fre¬ 
quency  is  shown  in  Fig.  4.  Variations  in  the  Young’s  modulus  of 
the  bi-polar  plates  and  the  GDEs  both  have  the  same  impact  on 
the  natural  frequency  of  the  material,  with  increases  resulting  in 
a  linear  increase  to  the  natural  frequency.  With  a  25%  increase  in 
the  Young’s  modulus,  the  natural  frequency  increases  by  approxi¬ 
mately  1 2%.  On  the  other  hand,  the  natural  frequency  is  not  affected 
by  increasing  the  Young’s  modulus  of  the  membrane  since  its  thick¬ 
ness  is  significantly  smaller  than  the  other  two  layers. 

4.3.  Influence  of  material  density 

Fig.  5  illustrates  the  effect  of  independently  increasing  the  den¬ 
sity  of  the  bi-polar  plates,  GDEs  and  membrane  on  the  lowest 
natural  frequency.  The  combination  of  the  order  of  magnitude 
of  density  and  thickness  will  be  the  determinant  factors  in  this 
behaviour.  Since  the  densities  of  the  bi-polar  plate  and  mem¬ 
brane  are  almost  4.5  times  greater  than  the  GDE,  and  the  thickness 


Fig.  5.  Effect  of  independently  varying  the  component  density  on  first  natural  fre¬ 
quency. 


of  the  GDE  and  membrane  are  much  less  than  the  thickness  of 
the  bipolar  plate,  increasing  the  density  for  the  GDEs  and  mem¬ 
brane  do  not  result  in  any  significant  changes  to  the  lowest 
natural  frequency.  On  the  other  hand,  increasing  the  density  of 
the  bi-polar  plates  results  in  a  strong  linear  dependence  of  the 
lowest  natural  frequency.  Increasing  the  density  by  25%  results 
in  a  decrease  of  the  lowest  natural  frequency  by  approximately 
11%.  The  effects  of  this  parametric  study  have  also  been  per¬ 
formed  on  the  2nd  and  3rd  modes,  which  resulted  in  similar 
trends. 


5.  Conclusions 

A  three-dimensional  finite  element  model  for  the  PEMFC  mod¬ 
elled  as  a  laminated  composite  structure  using  Mindlin’s  plate 
theory  is  presented.  The  model  is  employed  to  characterize  the 
vibration  behaviour  of  the  PEMFC  modelled  with  the  membrane, 
GDEs,  and  bi-polar  plates  considered  as  composite  material  plies. 
A  base  case  simulation  is  performed,  and  the  mode  shapes  for  the 
first  three  modes  are  presented. 

A  parametric  study  is  performed  to  investigate  the  effect  of  inde¬ 
pendently  increasing  the  thickness,  Young’s  modulus,  and  density 
of  the  bi-polar  plates,  GDEs  and  membrane  on  the  lowest  natu¬ 
ral  frequency.  It  is  found  that  increases  in  the  GDE  and  membrane 
thickness  have  little  effect,  while  increases  in  the  bi-polar  plate 
thickness  has  the  most  dominant  effect  on  the  natural  frequency  of 
the  entire  study.  To  a  lesser  degree,  increasing  the  Young’s  mod¬ 
ulus  of  the  bi-polar  plates  and  GDEs  and  increasing  the  density 
of  the  bi-polar  plate  also  have  a  significant  effect  on  the  lowest 
natural  frequency.  Increasing  the  Young’s  modulus  of  the  bi-polar 
plates  and  GDEs  result  in  a  linear  increase  in  the  natural  frequency, 
while  increasing  the  density  of  the  bi-polar  plates  results  in  a  linear 
decrease  in  the  natural  frequency.  These  results  provide  insight  into 
how  the  natural  frequencies  of  the  PEMFC  may  be  tuned  to  avoid 
high  amplitude  vibrations  by  modifying  the  material  and  geometric 
properties  of  individual  components. 

The  magnitude  of  the  natural  frequencies  for  this  single  PEMFC 
have  been  found  to  be  out  of  the  range  of  frequencies  (1-25  Hz) 
encountered  in  common  transportation  and  stationary  applica¬ 
tions  [1-3].  Future  work  includes  the  representation  of  the  PEMFC 
as  separate  layers  pressed  together,  so  that  the  interfacial  effects 
between  layers  can  be  investigated.  In  this  investigation,  one  sin¬ 
gle  cell  is  investigated  to  provide  insight  into  the  dominance  of  the 
component  layers.  However,  practical  applications  require  fuel  cell 
stacks  to  provide  realistic  power  levels.  Since  free  vibrations  that 
coincide  with  the  natural  frequency  will  result  in  high  amplitude 
resonance  vibrations,  further  investigations  are  required  to  deter¬ 
mine  whether  the  fuel  cell  stack  will  experience  failure  modes  due 
to  these  application-based  vibrations.  Furthermore,  in  these  practi¬ 
cal  applications,  fuel  cells  will  be  assembled  under  compression,  so 
our  future  investigations  will  include  multiple  PEMFC  cells  assem¬ 
bled  together  to  form  a  fuel  cell  stack,  which  will  provide  the  proper 
context  in  which  to  investigate  vibration-resistant  material  prop¬ 
erties  and  stack  geometry.  Furthermore,  the  absence  of  flow  fields 
in  the  bi-polar  plate  structure  may  also  affect  the  results  presented 
here;  therefore  the  flow  field  will  also  be  included  in  our  future 
work. 
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